V 


/<£•%/ at stJ / 


/ /A- 


■ — ✓ 


99098 

91 ? 


FINAL REPORT 


NASA GRANT NAG3-368 
1 October 1982 - 30 September 1986 


The Temperature Dependence of Inelastic Light Scattering 

from Small Particles 

for Use in Combustion Diagnostic Instrumentation 


prepared by 

Stanley D. Cloud 
Associate Professor of Physics 
University of Nevada, Las Vegas 
Las Vegas, Nevada 89154 


submitted to 


The National Aeronautics & Space Administration 
Lewis Research Center 
Cleveland, Ohio 


October 9, 1987 


(NASA-CB- 180399 ) THE TEMEEBAICBE DEPENDENCE 
Cf INELASTIC LIGHT SC A11EB ING EEC* SHALL 
PABTICLES fOfi DSE IN CCMBIS11CN DIAGNOSTIC 
1NSTB OHENTAIICN final Bepcrt, 1 Cct. 1982 
3C Sep* 1586 (Nevada Cniv.) 27 p Avail: 


N87-28634 


Unclas 

0CS90S8 



INTRODUCTION 


I . 


The work described In this report falls Into two parts: (1) 

a computer calculation of the expected angular distribution of 
coherent anti-Stokes Raman scattering (CARS) from micrometer size 
polystyrene spheres based on a Mie-type model, and (2) a pilot 
experiment to test the feasibility of measuring CARS angular 
distributions from micrometer size polystyrene spheres by simply 
suspending them in water. The computer calculations predict very 
interesting structure in the angular distributions that depends 
strongly on the size and relative refractive index of the 
spheres. The expected sharp structure makes an experiment 
particularly important not only as a test of the theoretical 
model's details but also as preliminary to a variety of possible 
applications . 

This work grew out of a project to find means for increasing 
the sensitivity or applicability of CARS as a tool for combustion 
diagnostic instrumentation. A great deal of work has been done 
in the 31 DO 1 1 cat ion of CARS to nnmhiistlnn stiid 1 1 , 2 vet man\r 

* r — * • ^ — - Jl 

problems remain to be solved. The theoretical work of Cooney and 
Gross 3 in 1982 showed that CARS from small particles might be 
enhanced by an order of magnitude over the intensity of 
scattering from the bulk material and by an even greater amount 
over the intensity of scattering from gas molecules. This 
suggested the possibility of seeding a flame with particles of 
suitable material and using CARS from the particles as a probe of 
local gas density and local pressure. 

Since only Cooney's and Gross's computation of CARS from 
benzene droplets had been made and since no experiments had been 
performed at all, it was essential to develop a basis of 
scientific information about particle CARS before passing 
judgment on the technique as a combustion diagnostics 
instrumentation tool. The first step was to extend the 
computational capability to include any material and to produce 
complete angular information about CARS from a particle of any 
size within the limits of the Mie-type calculation. The second 
step was to select a convenient particle material and attempt to 
measure the CARS angular distribution for comparison with the 
calculational model. 


II. THEORETICAL CALCULATIONS 

The first reported computer calculations of CARS from 
micrometer size particles were by Cooney and Gross 3 . Due to 
limitations of computer facilities, they only considered CARS 
from benzene liquid droplets, and they only did computations for 
thirteen different droplet sizes. Earlier experiments* with 
fluorescence from dye-doped latex spheres had shown sharp "size 
resonances" that might appear in CARS also. In order to indicate 
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the presence of such resonances, which are less than one percent 
In width, calculations needed to be made for a range of closely 
spaced particle sizes. The Cray-lS supercomputer at the Lewis 
Research Center made such detailed calculations possible. 

Details of the Mie-type model used for the computation of 
CARS angular distributions from particles are contained in the 
Ph.D. dissertation of Gross . 5 The essential ideas are the 
following. In CARS, three incident photons interact 
simultaneously with a molecule to produce a fourth photon with a 
different energy. Typically, two of the incident photons are 
frequency-doubled Nd : YAG laser photons with a wavelength of 
532 nm, while the third incident photon comes from a tunable dye 
laser driven also by the Nd : YAG laser. The dye laser wavelength 
is chosen in accordance with the Raman resonance one desires to 
use in a particular material. If the energy of the Nd : YAG laser 
photons is and the energy of the dye laser photon is E2# then 
the fourth photon resulting from CARS will have energy ZE^^-Ez- 

From the point of view of electric fields, in the CARS 
process described above, two of the incident photons induce 
fields in the particle that oscillate at frequency fj while the 
third induces a field oscillating at frequency f 0 . The presence 
of these fields simultaneously superimposed in the material 
produces a new electromagnetic oscillation at a frequency f 3 
equal to The new frequency results because the material 

has a third-order, nonlinear susceptibility. 

In order to predict the intensity and directional 
distribution of the new CARS produced light theoretically, the 
oscillation at frequency f 3 =2f2~fi is introduced as a source into 
Maxwell's equations. Maxwell's equations combine to form a wave 
equation with a source term at frequency f 3 whose intensity is a 
function of position within the particle. The wave equation can 
be solved numerically by the Green function technique and 
multipole expansion, subject to boundary conditions at the 
surface of the sphere. A mathematical outline of the procedure 
is given in appendix A. 

The computation routine involves several steps. First, the 
field distributions at frequencies f, and f2 are found inside the 
sphere in terms of the amplitudes ana polarizations of the 
incident photon beams. To do this, the fields are expressed in 
terms of Ricat t i-Bessel functions and associated Legendre 
functions. Then the distribution of sources at f 3 is calculated 
and integrated with the vector multipoles over the volume of the 
particle. The vector multipoles are computed from spherical 
Bessel functions and irreducible tensors, which are expressed in 
terms of spherical harmonics by means of the Clebsch-Gordan 
coefficients. Next, the integrals over the volume are combined 
with the boundary conditions at the surface of the sphere to 
yield coefficients in a multipole expansion of the scattered 
field. The scattered field is computed angle by angle as the 
partial sum of the expansion, including enough terms to give 
sufficient precision in the final result. 
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The program Is exhibited In appendix B. It Is very simple to 
use, because very little Input data are required. The size 
parameter of the particle Is Input at line 3300, and the number 
of terms to be Included In the multipole expansion at line 3400. 
The three wavelengths and the corresponding refractive indices 
are entered directly into the common block as REFRAC(3) and 
WAVE1(3). The entered wavelengths are those outside the 
particle. The program computes the wavenumbers inside the 
particle at line 4000. 

In operation, the program begins calculating the incident 
fields inside the particle using subroutines MIEC0 and INFLD. 
MIEC0 is called at line 4200. The special functions are called, 
and then at line 6400, the integration of the new field over the 
particle volume is begun. DPR0D (line 8700) is the new field at 
f 3 . The resulting integrals are CLE and CLM, and they are the 
coefficients to be used in the multipole sums begun in line 
14200. Output of the program is an angular distribution of the 
logarithm of the CARS intensity in relative units, printed at 
line 18600. FIT is the scattering angle in degrees measured from 
the direction of the incident beams. DCS is the log intensity. 

The subroutines used are: 

MIECO calculates the Mie-type coefficients for use 
in finding the internal fields in INFLD 

RBES computes the Ricatt i-Bessel functions; 

RHANK computes the Ricatt i-Hankel functions; 

INFLD finds the internal fields; 

LEG calculates the Legendre functions; 

YLM computes the spherical harmonics; 

TLM computes the irreducible tensors or vector 
spherical harmonics; 

CLG calculates the Clebsch-Gordan coefficients. 

Note that wherever possible the program is vectorized to run 
with optimal speed on a Cray-1 computer. For LMAX = 20, which is 
sufficient for size parameters up to about 8 , the program runs in 
approximately 40 sec of CPU time. Run time increases rapidly 
with increasing LMAX. With LMAX = 95, the program took about 3.3 
hours . 

CARS angular distributions were computed on the NASA Lewis 
Research Center Cray-lS computer, transmitted via Telenet to an 
IBM-PC , and stored on floppy disks. The angular distributions 
were then sent by telephone line from the IBM-PC to the 
University of Nevada, Las Vegas Computing Center Cyber-172 
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computer where they were archived and plotted using the National 
Center for Atmospheric Research graphics package called NCAR . 

The first calculations were a set of angular distributions 
for 29 different size parameters ranging from 5.39 to 6.79 for 
a polystyrene sphere Immersed in air. These results are 
presented In figures 1, 2, and 3. The vertical scale is the 
logarithm of the scattered intensity on a relative scale. 

Polystyrene was chosen as the material in anticipation of an 
experiment. Polystyrene is a polymer of styrene, whose structure 
includes a benzene ring. So polystyrene has Raman scattering 
properties very similar to benzene; in particular, it exhibits a 
strong, isolated Raman resonance near 1000 cm - !. 6 Polystyrene 
spheres carefully manufactured in a range of sizes are available 
commercially as "uniform latex". 

Figures 1 through 3 show angular distributions with 
considerable structure that changes rapidly at certain scattering 
angles as the particle size changes. The change of structure 
with size is a property that may be exploitable for applications. 
This will be discussed further in the section entitled 
"Discussion". However, the figures show that the change of 
angular distribution with size is very complicated, so great 
detail will be required for applications. 

Figure 4 shows the variation of CARS angular distributions as 
the refractive index changes. The refractive index is, of 
course, the relative refractive index, so a change can be the 
result of something happening within the particle or in the 
surrounding medium. The figure shows that CARS scattering 
changes relatively slowly with refractive index. Figure 5 shows 
that there is practically no change in CARS angular distribution 
when the ambient air pressure is increased from zero to 90 
atmospheres . 

For the purpose of doing an experiment to compare actual CARS 
with the computed angular distributions, one wants a scatterer 
that will produce structured distributions that are not sensitive 
to experimental conditions. For example, in a sample of 
polystyrene latex there might be a variation of 1* or 2 % in the 
sizes of the particles. If the scattering is sensitive to such 
small changes, then the particle distribution will smear the 
angular distributions. Providing the size distribution is known, 
it is possible to integrate over the range of angular 
distributions to obtain a smeared theoretical distribution. 
However, one should expect violently changing conditions in the 
vicinity of the high-powered laser beam required for CARS, and 
both the nature and the magnitude of the changes may be 
difficult to determine precisely. For example, due to heating or 
ablation, the actual size of the particles might change 
appreciably. Furthermore, the surrounding medium could expand or 
even vaporize if it is a liquid. In general, when the relative 
refractive index is reduced, the amount of structure and the 
size-variation of the structure in CARS angular distributions is 
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reduced. This Is seen In figures 6 and 7, which present CARS In 
water for a range of size parameters and for a range of 
refractive Indices. It suggests that for a preliminary 
experiment, It would be desirable to suspend the particles In a 
medium, such as water, whose refractive Index Is different from 
that of polystyrene (approximately 1.6), so that some structure 
will be present. But the surrounding refractive index should not 
be too different. A relative refractive Index of about 1.2, as 
shown In the figures, should give good results. 

These computations were presented in a conference paper. 7 


III. EXPERIMENT 

Preliminary experiments to observe CARS from polystyrene 
latex spheres were carried out using the Nd:YAG/Dye laser CARS 
system at the NASA Lewis Research Center. The plan included 
generating CARS from latex suspended in water held in a cell 
within a Wyatt Technology "Dawn Model B" fifteen-detector light 
scattering system. 8 The fifteen detectors were intended to 
detect the scattered light at fifteen angles simultaneously. 
Output from the detectors would be digitized and stored in a 
personal computer for later analysis. Interference filters 9 to 
transmit the CARS generated light but block directly scattered 
incident light were fitted in front of the fifteen detectors. 
Although the experiment did not proceed to the point of producing 
CARS within the "Dawn" system, the "Dawn" system was interfaced 
to the computer and tested satisfactorily as a detector system. 
Furthermore, the interference filters were found to discriminate 
very well against light from the laser beams, transmitted the 
CARS wavelengths as expected, and were extensively used in the 
pilot CARS study to be described. 

Before attempting to observe CARS from polystyrene particles, 
an effort was made to see CARS from bulk liquid benzene. A 
quartz cylindrical cell filled with benzene was positioned with 
its flat end faces perpendicular to the approximately colinear 
incident beams. With the dye laser wavelength set at about 562nm 
and its intensity roughly one-half the intensity of the Nd.-YAG 
laser beam at 532nm, a spot of blue light could be seen slightly 
off axis when the light exiting the cell was allowed to fall on a 
white card. It was necessary to view the card through the 
interference filters mentioned above in order to suppress the 
incident laser beams and also the yellow and red Raman shifted 
laser beams. The blue light was directed into a spectrometer via 
a prism, and the wavelength was measured to be about 505nm in 
agreement with the f 3 = 2 f 2 - fi expectation. 

To explore the possibility of detecting CARS from polystyrene 
microspheres, a very dilute suspension of 1.091 micrometer 
latex particleslO in water was prepared. The concentration of 
polystyrene in the cell was about 5 x 10-7, which would give 
approximately 5 microspheres in the interaction region where the 
incident laser beams were focused. The dimensions of the 
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interaction region were estimated to be 1mm of length by 0.1mm 
diameter. Using the Interference filters to block scattered 
incident laser light, it was possible to observe visually blue 
light emerging from the region of beam focus. A fiber bundle was 
used to collect the blue light and direct it through a lens that 
focused it on the slit of a double spectrometer. To confirm that 
it was CARS from the polystyrene spheres, the light passing 
through the spectrometer was analyzed by an optical multichannel 
analyzer system gated synchronously with laser pulses. As 
expected, a wavelength peak was found at about 505nm. However, 
the 505nm light remained present as the dye laser wavelength was 
scanned over a lOnm or more range, whereas a peak with a width of 
perhaps 20cm -1 (see ref. 6) or about 0.6nm was expected. This 
"continuous" background was attributed to weak CARS from a 
distant Raman resonance in water. To detect the presence of 
polystyrene, then, it was necessary to search for a weak, narrow 
peak superimposed on the water background in a plot of 505nm 
light intensity versus dye laser wavelength. 

Data collection was quite difficult because it was necessary 
to keep the laser beam intensities down so as to minimize damage 
to the cell window and also prevent boiling of the water. This 
meant that the signal was weak and difficult to collect 
efficiently with the optical fiber bundle. Integration periods 
of several minutes were required, during which the laser output 
power could drift, so laser power was monitored and regulated 
manually. Furthermore, whenever window damage was detected by the 
presence of photoacoustic pulses, the cell had to be moved 
slightly and the fiber bundle had to be repositioned. Average 
laser power was typically 13mW at 532nm and 6mw at 562nm with a 
10 hertz repetition rate. 

Late on the next to last day of scheduled time, just such 
a narrow peak was found in the day's last run. The data are 
shown in figure 8. Both peak position and width are as expected. 
Unfortunately, due to the difficulties mentioned above, it was 
not possible to confirm the presence of the peak. 


DISCUSSION 

The theoretical results indicate that CARS from particles is 
rich in structure. Whether the structure can be exploited 
efficiently for instrumentation purposes remains to be 
determined, and much more developmental work is needed. In 
particular, the fundamental question of the applicability of the 
Mie-type model to droplets of benzene or to polystyrene spheres 
has not yet been answered experimentally. The particle's surface 
is included in the model as an interface at which boundary 
conditions must match. The matching of boundary conditions 
requires the addition of a field in addition to that directly 
attributable to the sources at f 3 , and that fle ld results from 
the surface. But the model does not take into account that the 
CARS interaction might be different at the surface or that 
surface currents might be produced. Especially with such small 
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particles with large surface to volume ratios under the influence 
of very intense incident fields, surface effects may be 
significant. The Yale group has been successful in observing 
nonlinear scattering from droplets recently 11 , but they are 
working with larger droplets, of about 30 micrometers, rather 
than the approximately 1 micrometer particles considered here. 

Some modifications in the apparatus should be made before 
further experiments are attempted. If possible, laser power 
stability should be improved. In any case, the scattered light 
collection system should be made more efficient. Perhaps CARS 
signal strength can be increased by using a dense aerosol of 
particles Instead of a dilute suspension. 

Although the experiment described here was conceived more 
optimistically than it turned out, it was nevertheless intended 
only as a pilot study to determine the feasibility of detecting 
CARS from small particles and to discover what the main barriers 
to success are likely to be. It was successful in those 
respects: CARS was observed, if only briefly from polystyrene 

particles; and a clear idea of the difficulties has emerged. 
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APPENDIX A 


Mathematical Outline of the Procedure 
for Computing CARS from 
Spherical Dielectric Particles 


I. The outgoing CARS power per unit solid angle Is computed from 
the electromagnetic field amplitude, H: 

d p r c 

— ( 0 3 <p) = “I. r2 [ Hg »1?3 * ] » where 

dXi. atf 33 


II. H Is computed as an expansion of irreducible tensors: 

e lk3r 

^(-i) 1 + 1 ^CM(l,m)U x T llm ) + CE(l,m)T llm 


H 3 = 


7 * 


A 


T. . 

- urn 


LA 


C( 111 !A ) Yi m w (4 , 

* ^ ' l ■ 


III. and the coefficients are: 


(xh) ' ( yj ) - n 3 (yj ) ' (xh) 


CE( 1 ,m) 


n 3 (xh) ' (yj ) - (yj ) ' (xh) 


with 




(E) 


j 1 ( n 3k3 r ^ T llm( - 0 - , <J> ) 

- ---- V x ^lm ln < M )' 

n 3 k 3 


and 




n is the refractive index, k the wave number. h and j are 
the spherical Hankel and Bessel functions. x is k 3 a, y is 
n 3 X, and a is the particle radius. 
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APPENDIX B 


THE CRAY-1 PROGRAM 


"SCARS" 



0000800 PROGRAM SCARS 

0000900 COMMON/PROPS/P I . PMB » R A » REFRAC < 3 ) » WAVE 1 < 3 ) * WAVE2 (3) .LMAX 

000 1 000 COMMON/COORDS/R » COST . S I NT » COSP » S I NF » P 

0001100 COMMON/F1ELDS/ER<2)»£T(2) »EP(2) 

0001200 C0MM0N/RAD2/RB < 100 ) » DRB ( 100 ) 

0001300 C0MM0N/RAD3/RH ( 1 00 ) > DRH <100) 

0001400 C0MM0N/RAD4/RB2< 100) .DRB2< 100) 

0001500 COMMON/THETA/PIE ( 100) * TAU( 100) 

0001600 COMMON/COEFS 1 / AM I E < 1 00 » 2 ) » BMI E < 1 00 » 2 ) 

0001700 CQMM0N/C0EFS2/CMIE ( 100 . 2) * OMIE < 100* 2 ) 

0001800 COHMON/HAR/Y < 100 » 100) 

0001900 COMMON/CGC/C < 3 » 3 » 100 » 100) 

0002000 COMMON/TEN/TB < 3.3.100.100) 

0002100 COMPLEX ER.ET.EP 

0002200 COMPLEX TB 

0002300 COMPLEX RH * DRH 

0002400 COMPLEX AMIE.BMIE 

0002500 COMPLEX CNIE.DHIE 

0002600 COMPLEX Y 

0002700 COMPLEX SMI (3) »SM2<3) »SE(3) >SDB( 3) 

0002800 COMPLEX DM< 100) .DE<100) 

0002900 COMPLEX CLM( 100. 100) .CLE< 100.100) 

0003000 COMPLEX AM ( 3 » 100 » 100 ) * AE ( 3 , 100 . 100 ) » APM ( 100 » 100 ) . APE ( 1 00 . 100 ) 

0003100 COMPLEX DFROD, IMAGI >CIL 

0003200 DIMENSION VL< 100) .VPL< 100) 

0003300 SIZEP=70i0 

0003400 LMAX=95 

0003500 LLMAX=LMAX-1 

0003600 CALL CLG 

0003700 DO 1 1=1.3 

0003800 WAVE 1 ( I ) =2 . 0*P I /WAVE 1(1) 

0003900 C ABOVE CONVERTS INPUT WAVELENGTHS TO WAVENUMBERS 

0004000 1 WAVE2 < I ) =UAVE 1(1) JREFRAC ( I ) 

0004100 RA=SIZEP/WAVE1 ( 1 ) 

0004200 CALL MIECO(l) 

0004300 CALL MIEC0<2) 

0004400 C SAVE CARS "INSIDE* RICATTI BES. FUNCS. AS RB2 S DRB2 
0004500 ARG=WAVE2(3)*RA 

0004600 CALL RBES(ARG) 

0004700 DO 2 L=t»LMAX 

0004800 RB2(L)=RB(L) 

0004900 DRB2(L)=DRB(L) 

0005000 DL=1 .0/FLOAT (L+L+l ) 

0005100 FL=FLOAT (L) 

0005200 VL(L)=SQRT(FL*DL) 

0005300 VPL(L)=<FL+1.)*DL 

0005400 VPL(L)=SQRT(VPL(D) 

0005500 2 CONTINUE 

i.) 0005600 C CREATE CARS "OUTSIDE" RIC. HANK. FUNCS. 

0005700 ARG=UAVE1(3)*RA 

0005800 CALL RBES(ARG) 

O 0005900 CALL RHANK(ARG) 

i 0006000 3000 FORMAT! 18H BEGIN INTEGRATION) 

0006100 C 

O 0006200 CBEGIN INTEGRATION TO CREATE CLM(L.H) I CLE(L.M) 

0006300 C 




u 



0006400 

DO 10 11=21*40 

0006500 

3100 FORMAT (5H 11= *15) 

0006600 

X=0 . 05IFL0AT ( I 1 ) - 1 . 025 

AAAA 7AA 

XS=X*X 

0006800 

DO 10 12=21*40 

0006900 

3200 F0RMAT(5H 12= *15) 

0007000 

CY=0.05*FLOAT<I2)-1.025 

0007100 

YS=CY*CY 

0007200 

DO 10 13=1*40 

0007300 

Z=0.05*FL0AT<I3)-1.025 

0007400 

ZS=Z*Z 

0007500 

RS=XS+YS+ZS 

0007600 

IFtRS.GT. 0.998) GOTO 10 

0007700 

IFIRS.LT.l.E-4) GOTO 10 

0007800 

R=SGRT(RS) 

0007900 

COST=Z/R 

0008000 

RO=SQRT(XS+YS) 

0008100 

SINT=RO/R 

0008200 

COSP=X/RO 

0008300 

SINP=CY/RO 

0008400 

P=ASIN(SINP) 

0008500 

CALL INFLD(l) 

0003600 

CALL INFLDI2) 


- • :■*» / 3 r* ■' - 

’t >*_. »w 

OF POOR QUALITY 


0008700 BPR0D=ER( 1 )*C0NJG(£R(2) )+ET ( 1 )*C0NJG(ET (2) )+EP( 1 )*C0NJG(EP(2) ) 

0003800 C GENERATE A(L*M)S 

0008900 ARG=WAVE2<3)*R*RA 


0009000 CALL RBES(ARG) 

0009100 DO 3 LD=1 *LMAX 

0009200 3 RB(LD)=RB(LD)/ARG 

0009300 C ABOVE MADE SPHER. BESS. FUNC. FROM RICATTI B. F. 

0009400 RBO=S I N ( ARG ) / ARG 

0009500 CALL YLM 


0009600 CALL TLM 

0009700 C AE S AM HERE ARE MULTIPOLE VECTORS* NOT COEFFS. 

0009800 DO 5 IC=1»3 

0009900 AM(IC*1»1)=TB<IC,2*1,1)*RB<1) 

0010000 AN(IC.1*2>=TB<IC»2.1*2)*RB<1) 

0010100 AE(IC,l*l)=-VL(imB<IC»3»l*l>*RB<2>+VPL<lWB(IC*l*l*l)*RB0 

0010200 5 AE(ICfl>2)=-VL(imB(IC»3fl>2)*RB(2)FVPL(l)*TB(ICfl»lf2)*RB0 

0010300 DO 9 L=2»LLMAX 

0010400 MT0P=LF1 

0010500 DO 8 M=1 »MTOP 

0010600 AM(1,L*M)=TB(1,2»L*M)*RB<L) 

0010700 rtM(2*L*M)=TB(2*2*L*M)4RB(L) 


0010300 AM<3*L*H)=TB<3*2*L*H)*RB<L) 

0010900 AE( 1 *L*M)=-VL<L >4TB( 1 *3*L*H)IRB(L+1)+VPL(L)*TB( 1 * 1 *L*M)#RB(L-1 ) 

0011000 AE<2»L»M>=-VL<L>*TB<2*3»L*M)*RB<LH>+VPL(L>*TB<2*1*L*M)*RB(L-1> 

0011100 AE<3*L»M>=-VL<L)*TB<3i3*L*M)*RB(L+lHVPL<L>*TB<3*l*L»N)*RB<L-l> 

0011200 8 CONTINUE 

0011300 9 CONTINUE 


0011400 C BEGIN NEW L LOOP STARTING WITH L=1 

0011500 DO 19 L=1 *LLMAX 

0011600 MT0P=L+1 

0011700 DO 18 M=1 iMTOP 

0011800 AM(l»L»M)=CONJG(AM(l»L*H) ) 

0011900 AH(2*L*M)=C0NJG(AH(2*L*H) ) 

0012000 AM(3>L*M)=C0NJG(AM(3tL>M) ) 

0012100 AE(lfL*M)=CONJG(AE(l»LiM)) 

0012200 AE(2fL?M)*C0NJG(AE(2»l»H)) 

0012300 AE(3»ltH)-C0NJG(AE(3'lfM)) 

0012400 18 CONTINUE 

0012500 DO 28 H-lfHTOP 

0012400 MHa»N)-(M(ltl»N)tCMl)4M(2tLfK)8n(i)4AH(l>LiN)tEP(l»tDPftOD 

0012700 APE<l»H>*CAEU»liH>tCR<mAE<2»lfH>tET<mAE<3>l»H>*EPU>>W*08 

0012800 ' 28 CONTINUE — . 

0012900 C CLN t CUE AflE. THE INTEGRALS CALLED BIGHA(M) I SIGMAtE) IN NOTES v 


W” 


0013000 DO 38 M=1»MT0F' 

0013100 CLM<L»M)=CLM<L>M)+APM<L>M) 

0013200 CLE ( L » M ) =CLE ( L » M ) TAPE ( L f N ) 

AAITXAA IQ rnWTTMIIC 

VVtUWVV W W w w i I 4I1UL 

0013400 19 CONTINUE 

0013500 10 CONTINUE 

0013600 C CLE ( L > M ) =C0N JG < CLE ( L > M ) ) 

0013700 C CLM(L»M)=CONJG(CLM(L»M) ) 

0013800 1500 F0RMAT(E15.7»3XfE15*7) 

0013900 C 

0014000 C BEGIN S SUM FOR EACH OBSERVATION ANGLE 
0014100 C 

0014200 IMAGI=<0.0»1.0> 

0014300 SINP=1.0 

0014400 C0SP=0.0 

0014500 P=PI/2. 

0014600 WRITE (12» 2000) 

0014700 DO 30 L=1 >LLMAX 

0014800 DM ( L ) =DRH ( L ) *RB2 ( L ) -REFRAC ( 3 ) *DRB2 ( L ) *RH ( L ) 

0014900 30 DE(L)=REFRAC<3>*DRH(L>*RB2<L>-RH<L>*DRB2(L> 

0015000 DO 70 IT=1»179 

0015100 FIT =FL0AT( IT) 

0015200 T=FIT*PI/180. 

0015300 SINT=S1N(T) 

0015400 C0ST=C0S(T) 

0015500 DO 40 IC=1,3 

0015600 SDB(IC)=(0.0>0.0) 

0015700 40 CONTINUE 

0015800 CIL=-IMAGI 

A A 4 rn A A r A I « \I I w 

VVIJ 7W LrMLU 1 1.11 

0016000 CALL TLM 

0016100 DO 60 L=1 »LLMAX 

0016200 DO 55 IC=1»3 

0016300 SM1UC) = <0. 0.0.0) 

0016400 SM2(IC)=(0.0»0.0) 

0016500 55 SE<IC>=(0. 0,0.0) 

0016600 MT0P=L+1 

0016700 DO 50 M=2»MT0P»2 

0016800 SM1{1)=SM1(1)+CLM(L»MXAIMAG(TB(1» 3»L»M) ) 

0016900 SMl(2)=SMl(2H-CLM(LfM)*AIMAG(TB(2»3»L»M)) 

0017000 SM1(3)=SM1(3)+CLM(L»M)*AIMAG<TB(3»3»L»M) ) 

0017100 SM2( 1 )=SM2( 1 )+CLMCL»M)4AIMAG(TB( 1 f 1 »L»M) ) 

0017200 SM2(2)=SM2(2)+CLM(L»M)*AIMAG(TB(2»1>L,M)> 

0017300 SM2(3)=SM2(3)+CLM(L»M)*AIMAG<TB(3fl>L»M)) 

0017400 SE( 1 )=SE( 1 )+CLE(L>M)*AIHAG(TB(l »2>L»M) ) 

0017500 SE(2)=SE(2)fCLE(L»M)*AIMAG(TB(2»2fL>M)) 

0017600 SE(3)=SE(3)+CLE<L»M)4AIMAG(TB(3»2»L»M) ) 
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0017700 50 CONTINUE 

0017800 CIL=-CIL*IMAGI 

0017900 DO 65 IC=1,3 

0018000 65 SDB(IC)=SDBaC)+CIL*((VL(L)*SMl(IC)+VPL(L)*SM2(IC))/DM(L)-IMAGI* 


0018100 CSE(IC)/DE(L>) 
0018200 60 CONTINUE 


0018300 1510 FORMAT (6(E15.7»3X) ) 

0018400 DCS=SDB ( 1 ) *C0N JG (SDB(l) ) +SDB ( 2 ) *C0N JG ( SDB ( 2 )HSDB < 3 > 1C0NJG ( SDB ( 3) ) 

0018500 DCS=AL0G10(DCS> 

0018600 WRITE ( 12 t 2100 )FIT» DCS 


0018700 70 CONTINUE 

0018800 1000 F0RHAT(F7.3» 1X» 13) 

0018900 1100 FORMAT (32H TYPE SIZEP(F7.3> SPACE LMAX(I3)> 

0019000 2000 FORMAT ( 16H(E15. 7» 3X*E15»7) ) 

0019100 2100 FORMAT (El5*7»3XfE15# 7) 

0019200 STOP f 

0019300 END i 

0019400 SUBROUTINE HIECO(IN) ' * 

OOlWftft COMMON/PROPS/PI >PMB»RA> REFRAC (3) »UAVF1 (3) »WAVE2(3) »LMAX 


0019600 
0019700 
0019800 
0019900 
0020000 
0020100 
0020200 
0020300 
0020400 
0020500 
0020600 
0020700 
0020800 
0020900 
0021000 
0021100 
0021200 
0021300 
0021400 
0021500 
0021600 
0021700 
0021800 
0021900 
0022000 
0022100 
0022200 
0022300 
0022400 
0022500 
0022600 
0022700 
0022800 
0022900 
0023000 
0023100 
0023200 
0023300 
0023400 
0023500 
0023600 
0023700 
0023800 
0023900 
0024000 
0024100 
0024200 
0024300 
0024400 
0024500 
0024600 
0024700 
0024800 
0024900 
0025000 
0025100 
0025200 
0025300 
0025400 
0025500 
0025600 
0025700 
0025800 
0025900 
0026000 
ArtOA* AA 


C0MMQN/RAB2/RB ( 1 00 ) * DRB ( 100 ) 
C0MM0N/RAD3/RH< 100) *DRH< 100) 
C0MM0N/RAD4/RB2 < 100 > * DRB2( 100 ) 


COnNuN/COEFSi/AnlEI 100*2) » BHIE< 100*2) 
C0MM0N/C0EFS2/CHIE( 100*2) *DMIE( 100*2) 
COMPLEX RH * DRH * AMIE . BMIE * CHI E * DMIE 
X=WAVE2(IN)*RA 
CALL RBES(X) 

DO 10 L=1*LMAX 


RB2(L)=RB(L) 

DRB2(L)=0RB(L) 


10 CONTINUE 

X=WAVE1(IN)*RA 
CALL RBES(X) 

CALL RHANK(X) 
XXX=REFRAC(IN)*PMB 
DO 20 L=1 *LMAX 
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AMIE<L*IN)=1.0/(PMB*RH(L)*DRB2<L)-REFRAC<IN)*DRH(L)*RB2(L)> 
DMIE(L*IN)=1.0/( PMBtDRHt L ) *RB2(L ) -REFRAC ( IN ) #RH ( L ) 4DRB2 ( L ) ) 
BMIE ( L * IN ) = ( REFRAC (IN) tRB < L ) 4DRB2 ( L ) -PMBtDRB ( L ) *RB2 ( L ) ) * 
CDMIE(L* IN) 

CMIE<L»IN)=XXX*(RB<L>*DRH(L)-RH<L)*DRB(L)) 

DMIE(L*IN)=CMIE<L*IN)»DMIE(L*IN) 

CMIE(L*IN>=-AMIE<L,IN)*CMIE(L*IN> 

AMIE(L*IN)=<REFRAC(IN)*DRB(L)*RB2<L>-PMB*RB<L)*DRB2(L>>* 

CAMIE(L.IN) 

20 CONTINUE 
RETURN 
END 


bUBKUUIlNt 

COMMON/PROPS/PI » PMB * RA * REFRAC(3) * WAVE1 ( 3 ) , UAUE2 ( 3 ) * LMAX 
C0MM0N/RAD2/RB < 100 > * DRB < 100 ) 

S=SIN(X) 

C=COS(X) 

RIC1=S/X-C 

IF(X.LT. 0.9)60 TO 25 
LBIG=LMAXf 10 
LBM1=LBIG-1 
RB(LBIG)=0. 1 
RB(LBM1)=0.2 
DO 10 L=2*LBM1 
LL=LBIG-L 
LLP1=LL+1 
LLP2=LL+2 
REALIE=2*LL+3 

RB(LL)=REALIE*RB(LLP1)/X-RB(LLP2> 

10 CONTINUE 


SCALE=RIC1/RB( 1 ) 

RB(1)=RIC1 

DRB(1)=S-RB(1)/X 


DO 20 L=2*LMAX 
RB(L)=SCALE*RB<L) 

LM1=L-1 

REALIE=L 

DRB(L) =RB (LM1) -REALIE*RB < L ) /X 
20 CONTINUE 
RETURN 

25 RB(1)=RIC1 

RB(2)=3.0*RB(1)/X-S 
DRB(1)=S-RB(1)/X 
DRB(2)=RB(l)-2.0tRB(2)/X 
FR0NT=X«3/15.0 
Z»0.5*X*X 
DO 40 L*3>LHAX 
REALIE-2BL+1 

rDnuT-rc>okiT*v /or ai tc 






0026200 


TERM=1 .0 

0026300 


SERIES=1.0 

b026400 


NUM=2*L+3 

AAOHAA 

V VbWWVV 


nn ta t-i . l 

i’W WV A — 4 r w 

0026600 


REALIE =NUM*I 

0026700 


TERM--TERM4Z/REALIE 

0026800 


SERIES=SERIES+TERM 

0026900 


NUM=NUM+2 

0027000 

30 

CONTINUE 

0027100 


RB<L)=FRONT*SERIES 

0027200 


LM1=L-1 

0027300 


REALIE=L 

0027400 


DRB(L) =RB ( LM1 ) -REALIE*RB ( L ) /X 

0027500 

40 

CONTINUE 

0027600 


RETURN 

0027700 


END 

0027800 


SUBROUTINE RHANK(X) 
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0027900 C0MM0N/RAD3/RH< 100) * DRH( 100) 

0028000 C0MM0N/RAD2/RB ( 1 00 ) » DRB ( 100 ) 

0028100 COMMON/PROPS/PI »PMB»RA»REFRAC<3) fWAVEl (3) »WAVE2(3) >LMAX 

0028200 DIMENSION RN( 100) t DRN( 100) 


0028300 COMPLEX RH»DRH 

0028400 S=SIN(X) 


0028500 C=C0S(X) 

0028600 RN( 1 )=-C/X-S 

0028700 RN(2)=3.0*RN(1)/X+C 

0028800 DRN ( 1 ) =-C-RN ( 1 ) /X 


0028900 DRN(2)=RN(1)-2.0*RN(2)/X 

0029000 RH( 1 )=CMPLX(RB( 1 ) >RN(1 ) ) 


0029200 DRH(l) =CMPLX < DRB < 1 ) > DRN < 1 ) ) 

0029300 DRH ( 2 ) =CMPLX < DRB < 2 ) > DRN ( 2 ) ) 


0029400 DO 10 L=3>LMAX 


0029500 LM1=L-1 

0029600 LM2=L-2 


0029700 REALIE=2*L-1 

0029800 RN ( L ) =REALIE*RN ( LM1 ) /X-RN ( LM2 ) 

0029900 REALIE =L 

0030000 DRN ( L ) =RN ( LM1 ) -REAL IE*RN ( L ) /X 

0030100 RH(L)=CMPLX(RB(L) »RN(L) ) 

0030200 DRH ( L ) =CMPLX ( DRB < L ) » DRN < L ) ) 

0030300 10 CONTINUE 

0030400 RETURN 

0030500 END 


0030600 SUBROUTINE INFLD(NW) 

0030700 COMMON/PROPS/P 1 1 PMB » RA t REFRAC < 3 ) t WAVE 1 ( 3 ) » WAVE2 ( 3 ) t LMAX 

0030800 COMMON/COORDS/R » COST » SI NT » COSP i SINP » P 

0030900 COMMON/FIELDS/ER < 2 ) > ET < 2 ) f EP ( 2 ) 

0031000 C0MM0N/RAD2/RB ( 1 00 ) » DRB ( 1 00 ) 

00311 00 COMMON/THETA/PIE < 100 ) , TAU < 100) 

0031200 C0MM0N/C0EFS2/CM I E < 1 00 , 2 ) » DMIE ( 1 00 > 2 > 

0031300 COMPLEX ER,ET,EP>CMIE»DMIE»C»IMAGI 

0031400 IMA6I=(0.0» 1 .0) 

0031500 RSIZE=WAVE2(NW)*R*RA 

0031600 CALL RBES<RSIZE) 

0031700 CALL LEG 

0031800 ER(NH)=(0.0»0.0) 

0031900 ET(NU)=(0.0f0*0) 

0032000 EP(NW)=(0.0»0.0) 

0032100 C=(l. 0»0.0) 

0032200 DO 10 L=1*LMAX 

0032300 A=2*L+1 

0032400 B®Lt(L+l) 

0032500 ER(NU)»ER(NU)iC«A*CHIE(LrNU)tRB<L)tPIE(L)tSINT 

0032600 A-A/B 

0032700 C»C*IMABI 



0032800 

0032900 

0033000 

0033100 

0033200 

0033300 

0033400 

0033500 

0033600 

0033700 

0033800 

0033900 

0034000 

0034100 

0034200 

0034300 

0034400 

0034500 

0034600 

0034700 

0034800 

0034900 

0035000 

0035100 

0035200 

0035300 

0035400 

0035500 

0035600 

0035700 

0035800 

0035900 

0036000 

0036100 

0036200 

0036300 

0036400 

0036500 

0036600 

0036700 

0036800 

0036900 

0037000 

0037100 

0037200 

0037300 

0037400 

0037500 

0037600 

0037700 

0037800 

0037900 

0038000 

0038100 

0038200 

0038300 

0038400 

0038500 

0038600 

0038700 

0038800 

0038900 

0039000 

0039100 

AAIO'JAA 


ET<NW)=ET<NW)+C*A*<DMIE<L*NU)*RB<L)*PIE<L>-IMAGI*CMIE(L*NW)* 

CDRB<L)*TAU(L)> 

EP(NW)=EP(NH)4C<!A*(!HAGIICMIE(L*NU)*DRB(L)»PIECL)-DMIE(L,NU)* 
CRB(L)*TAU(L) ) 

10 CONTINUE 

ER < NW ) =ER < NW ) *COSP/ ( RSI ZEM2 ) 

ET(NW)=ET(NW) *COSP/RSI ZE 
EP < NU ) =EP < NU ) *S INP/RSI ZE 
RETURN 
END 

SUBROUTINE LEG 

COMMON/PROPS/PI » PMB * RA * REFRAC ( 3 ) > WAVE 1 ( 3 ) » WAVE2 ( 3 ) » LMAX 
COMMON/COORDS/R i X r SINT » COSP * SINP * P 
COMMON /THETA/PIE( 100) »TAU(100) 

F'IE(1)=1.0 
PIE(2)=3.0*X 
TAU( 1 )=X 

TAU(2)=6.0*X*X-3.0 
DO 10 L=3iLMAX 
LM1=L-1 
LM2=L-2 
REAL1=L 
REAL2=LM1 
REAL3=2*L-1 

P I E ( L ) = ( RE AL 3*X*P I E < L M 1 ) -RE AL 1 *P I E < LM2 ) > /REAL2 
REAL2=L+1 

TAU(L)=REAL1*X*PIE<L)-REAL2*PIE(LM1) 

10 CONTINUE 
RETURN 
END 

SUBROUTINE YLM 
COMMON/HAR/Y ( 100 » 100) 

COMMON/COORDS/R * X * S » COSP » SINP > P 

COMMON/PROPS/P I * PMB * RA » REFRAC ( 3 ) » UA VE 1 < 3 ) * WAVE2 < 3 > » LMAX 
COMPLEX Y 
COMPLEX EIP(IOO) 

DIMENSION YY( 100>100) 

D I MENS I ON PTM < 1 00 > * SHP < 1 00 ) * CMP < 1 00 ) 

NMAX=LMAX+1 
MMAX=NMAX+1 
DO 1 1=1 * NMAX 
DO 1 J=1»MMAX 
1 YY(I?J)=(0.0*0«0) 

YY(1*2)=1. 

YY(2*2)=X 

YY < 2 * 3 ) =-S/SQRT ( 2 . ) 

YY(3»2)=(3.*X*X-l.)/2. 

YY(3f3)=-SQRT(1.5)*X*S 
Y Y ( 3 * 4 ) =SQRT ( 0 . 375 ) *S*S 
DO 100 L=3*LMAX 
FL=FLOAT(L) 

RL=FL-1. 

RLL=FL+RL 

N=L+1 

NL=N-1 

IL=NL-1 

YY < N » 2 ) = ( RLL<Xf YY ( NL * 2 ) -RLf YY ( IL * 2 ) ) /FL 
L2=L-2 

DO 200 M=1*L2 
NM=M+2 
NMM=NM-1 
FH=FLOAT(M> 

PM=FL4FM 

PMM=PM-1. 

DM=ri -CM 



0039400 

0039500 

0039600 

0039700 

0039800 

0039900 

0040000 

0040100 

0040200 

0040300 

0040400 

0040500 

0040600 

0040700 

0040800 

0040900 

0041000 

0041100 

0041200 

0041300 

0041400 

0041500 

0041600 

0041700 

0041800 

0041900 

0042000 

0042100 

0042200 

0042300 

0042400 

0042500 

0042600 

0042700 

0042800 

0042900 

0043000 

0043100 

0043200 

0043300 

0043400 

0043500 

0043600 

0043700 

0043800 

0043900 

0044000 

0044100 

0044200 

0044300 

0044400 

0044500 

0044600 

0044700 

0044800 

0044900 

0045000 

0045100 

0045200 

0045300 

0045400 

0045500 

0045600 

0045700 

0045800 

0045900 


Y Y < N f NM ) = ( SORT ( RH*RMM ) * YY ( IL > NM> -RLUS* YY (NLfNMM) ) /SORT ( PMIPMM ) 
200 CONTINUE 


L1=L-1 

r>n iaa y.i 4 . I 
uu n-uift 

NM=M+2 

NMM=NH-1 

PM=FL+FLOAT(M) 

PMM=PM-1. 

YY ( N f NM ) =-RLL*S*YY ( NL f NMM) /SORT ( PM4PMH ) 
300 CONTINUE 
100 CONTINUE 

FPI=0*5/SQRT(PI ) 

Y(1f2)=YY<1f2)*FPI 
DO 500 L=1 fLMAX 
NP=L+2 

RAD=FL0AT(L+L+1) 

FACT=FPI*SQRT(RAD> 


DO 800 NM=2fNP 
PTM(NM)=FL0AT(NM-2)*P 
800 CONTINUE 

DO 700 NM=2fNP 
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SMP(NM)=SIN(PTM(NM) > 
CMP(NM)=COS(PTM<NM) ) 

700 CONTINUE 

DO 900 NM=2fNF 

EIP(NM)=CMPLX(CMP(NM) fSMP(NM) ) 
900 CONTINUE 


DO 400 NM=2»NP 

Y(L+1fNM)=YY(LF1fNM>*EIP(NM)*FACT 

* a a nnvimmr 

muu uun 1 inuL 
500 CONTINUE 

DO 600 L=1 fLMAX 
600 Y ( L+ 1 f 1 ) =-C0N JG ( Y ( L+ 1 f 3 ) ) 

RETURN 

END 

SUBROUTINE TIM 

COMMON/PROPS/P I » PMB » RA , REFR AC ( 3 ) » WAVE 1 C 3 ) f WAVE2 (3) fLMAX 
COMMON/CGC/C < 3 f 3 f 100 f 100 ) 

COMMON/HAR/Y ( 100» 100) 

COMMON/TEN/TB < 3 f 3 f 100 » 100 > 

COHMON/COORDS/R f CT f ST f CP f SP f P 

COMPLEX TBf YfT(3f3f100f 100) f AIMfEIPfCEIPfSTPfCSTPfCTPfCCTP 
DO 2 NL=1 fLMAX 
MT0P=NLF1 
DO 3 NM=1 »MT0P 

T( 1 f 1 fNLfNM)=C( 1 f 1 jNLfNM) tY(NL»NM+2) 

T ( 1f2fNLfNM)=C( 1f2fNLfNM)4Y(NL+1 fNM+2) 

T ( 1 f 3 f NL f NM ) =C ( 1 f 3 f NL f NM ) >Y ( NL +2 f NM+2 ) 
T(2f1fNLfNM)=C(2f1fNLfNM)*Y(NL»NM+1) 
T<2f2fNLfNM)=C(2f2fNLfNM)*Y(NLHfNMH) 
T(2f3fNLfNM)=C<2f3fNLfNM)*Y(NL+2fNH+1) 

T(3f l fNLfNM)=C(3f 1 fNLfNM)IY(NLfNM) 

T(3f2fNLfNM)=C<3f2»NLfNM)*Y<NL+1fNM) 

T(3f3fNLfNM)=C(3f3fNL#NM)*Y(NL+2fNM) 

3 CONTINUE 
2 CONTINUE 
R2=0. 7071067 
AIM=(0.0f1,O) 

EIP=CMPLX(CPtSP)*R2 

CEIP=CONJG(EIP> 

STP=ST*EIP 

CSTP=CONJQ(STP) 

CTP»CT*EIP 
CCTP*C0N JO < CTP ) 

DO 5 NL»1»LMAX 
MT0P=NLfl 


0046000 

.0046100 

0046200 

0046300 

0046400 

0046500 

0046600 

0046700 

0046800 

0046900 

0047000 

0047100 

0047200 

0047300 

0047400 

0047500 

0047600 

0047700 

0047800 

0047900 

0048000 

0048100 

0048200 

0048300 

0048400 

0048500 

0048600 

0048700 

0048800 

0048900 

0049000 

0049100 

0049200 

0049300 

0049400 

0049500 

0049600 

0049700 

0049800 

0049900 

0050000 

0050100 

0050200 

0050300 

0050400 

0050500 

0050600 

0050700 

0050800 

0050900 

0051000 

0051100 

0051200 

0051300 

0051400 

EOF 

Q 


DO 6 NM=1 fMTOP 

TB< If 1 fNLfNM)=-STP*T(3f 1 fNLfNM)+CSTP<T ( 1 > 1 fNLfNM)+CT$T(2f 1 fNLfNM) 
TB(2f1fNLfNM)=-CTP*T<3f1fNLfNM>+CCTP*T<1f1fNLfNM)-ST*T<2f1fNLfNM> 


T£{ 3f 1 fNLfNM) =-AI n$(EIr*T(3F 1 fNLfNM) TCEIPST < If i>nL?Nn} ) 
TB<l»2iNL»NM)=-STP*T(3»2»NL»NN>fCSTP*T(l»2»NL»NMHCT*T<2»2»NLFNM) 


TB(2f2fNLfNM)=-CTP*T<3f2fNLfNM>+CCTP*T(1f2fNLfNM)-ST*T<2f2fNLfNM> 

TB<3f2fNLfNM)=-AIM*(EIP*T(3f2fNLfNM)+CEIP*T<1f2fNLfNM)) 

TB(1f3fNLfNM)=-STP*T<3f3fNLfNHHCSTP*T<1f3fNLfNM)+CT*T(2f3fNLfNM> 


TB(2f3»NL»NM)=-CTP»T(3»3.NLfNM)+CCTP»T(l,3*NL.N«)-ST*T(2»3fNL»NM) 

TB(3»3»NL»NM)=-AIM*(EIPfT(3»3»NLfNH)+CEIP*T(l»3»NL»NM)) 

6 CONTINUE 
5 CONTINUE 
RETURN 
END 

SUBROUTINE CLG 

COMMON/PROPS/PI t PMB f RA f REFRAC ( 3) f WAVE1 ( 3 ) f WAVE2 ( 3 ) f LNAX 
COMMON/CGC/C ( 3 f 3 f 1 00 f 100 ) 


DO 1 NL=1fLMAX 
L=NL 
NP=L+1 
FL=FLOAT(L) 
D1=FL+1 . 
D2=B1+FL+1 » 
D3=D2+1 . 
D4=2.*FL-1, 

DO 1 NM=1fNP 
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M=NM-1 

FM=FLOAT<M) 

FLM=FL-FM 
FPH=FLtFH 
FPM1=FPM+1 . 

FLM1=FLM+1. 

FPM2=FPM-1 . 

FLM2=FLM-1 . 

IF (L-M)2f2f3 

3 C ( 1 f 1 f NL f NH ) =SQRT ( FLM*FLH2/2 . /FL/D4 ) 

2 C ( 2 f 1 f NL f NM ) =SQRT ( FPM*FLM/FL/D4 ) 

C ( 3 f 1 f NLf NH ) =SQRT < FPM*FPM2/2 , /FL/D4 ) 

C ( 1 f 2 f NL f NM ) =SQRT ( FLM*FPM 1 /2 . /FL/D1 ) 

C ( 2 f 2 f NL f NM ) =FM/SQRT ( FL*D 1 ) 

C ( 3 f 2 f NL f NM ) =-SQRT ( FPM*FLM1 /2 . /FL/D 1 ) 

C ( 1 f 3 f NL f NM ) =SQRT ( FPMlf ( FPM1+1 , ) /D3/D2) 

C ( 2 f 3 f NL f NM ) =-SQRT ( FLM1*FPM1/D1/D3 ) 

C(3f3fNLfNM)=SQRT( (FLM1+1 . )*FLM1/D3/D2> 

1 CONTINUE 
RETURN 
END 

BLOCK DATA 

COMMON/PROPS/P I f PMB f RA f REFRAC ( 3 ) f NAVE 1 (3) fUAVE2<3) fLMAX 
DATA REFRAC/1 . 521 3 f 1 . 5170 f 1 .5255/ 

DATA WAVE 1/0.445f0.465f0. 426/ 

DATA PI/3.1415926536/ 

DATA PMB/l.O/ 

END 

/EOF 


o 


POLYSTYRENE IN AIR 
VARIATION OF CARS PATTERN 
AS SEE PARAMETER CHANGES 
FROM 5.39 TO 5.64 
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VARIATION OP CARS PATTERN 
AS SEE PARAMETER CHANGES 
FROM 0.84 TO 6.28 
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FIGURE 2 











:ll 


origjmal pag:: 

OF POOR QUALITY 






ORIGiJ'U -- ■> 

Of POOH QUALITY 



FIGURE 4 
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FIGURE 5 







POLYSTYRENE IN WATER 
VARIATION OF CARS PATTERN 
AS SEE PARAMETER CHANGES 
FROM 6.69 TO 7.59 
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FIGURE 6 











INTENSITY (arb units) 
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